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SUMMARY 

Calculations of Ice shapes and the resulting drag Increases are presented 
for experimental data on a NACA 0012 airfoil. They were made with a combina- 
tion of LEWICE and Interactive boundary-layer codes for a wide range of condi- 
tions which Include airspeed and temperature, the droplet size and liquid water 
content of the cloud, and the angle of attack of the airfoil. In all cases the 
calculated results account for the drag Increase due to Ice accretion and. In 

general, show good agreement with data. 

/ / 


1.0 INTRODUCTION 

The prediction of Ice shapes and the determination of their effect on 
lifting surfaces Is a problem of central Importance In aircraft design since 
ice accretion can adversely affect aerodynamic performance of aircraft compo- 
nents. In recent years, research has been undertaken to Improve understanding 
of the formation of the Ice, Its accretion and the consequences for aerodynamic 
performance. Reference 1 provides an overview of the analytical and experimen- 
tal Icing activities In progress and reference 2 a review of the progress on 
one element of the overall activity, namely the unprotected airfoil Icing prob- 
lem. The latter discusses the development and validation of computer codes 
which predict the buildup of Ice on unprotected airfoils and the resulting deg- 
radation of aerodynamic performance due to Ice accretion. 

These state-of-the-art reviews show that a two-dimensional Ice accretion 
code, LEWICE, developed In 1983 at the University of Dayton Research Institute 
(ref. 3) and later modified by Ruff (ref. 4), provides a basis for the determi- 
nation of the Ice buildup on the leading edge of airfoils, for both the glaze 
Ice formed at temperatures slightly below freezing and at relatively high liq- 
uid water contents and high flight speeds and for the rime Ice which occurs at 
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low temperatures, low liquid water contents, and low flight speeds. The three 
major elements of this code are: (1) an Invlscld panel code, (2) a water drop- 

let trajectory model, and (3) an energy balance equation proposed by Messinger 
(ref. 5). The Ice accretion Is computed on the airfoil leading edge as a func- 
tion of time with user specified time Intervals . At each given time, the 
Invlscld flowfleld Is determined from the panel code so that trajectory and 
heat transfer calculations can be performed. As the Ice accretion builds up, 
Its shape may become rugged especially In the case of glaze Ice which Is char- 
acterized by horns, and a rough, Irregular surface, and leads to higher aerody- 
namic losses unlike the rime ice. "Surface irregularities of the ice shape, see 
for example figure 1, can lead to multiple stagnation points with subsequent 
numerical difficulties, including a breakdown of the trajectory calculations 
which are necessary to determine the spatial distribution of water droplets. 

The automated smoothing procedure of reference 6 overcomes this difficulty by 
reducing the amplitude of the surface irregularities without loss of Important 
flow characteristics and usually allows the calculations to be performed for 
greater times than before, without the problems caused by multiple stagnation 
points. 

In addition to the need to predict Ice shapes accurately. It is important 
to determine the performance degradation of the airfoil due to icing. This can 
be achieved by two codes based on solutions of the Reynolds-averaged Navier- 
Stokes equations and their reduced forms. The Navier-Stokes method employs the 
ARC-2D code and has been developed for iced airfoils by Potapczuk (ref. 7) and 
the Interactive boundary-layer method of Cebeci (ref. 8) combines the solutions 
of the invlscid and viscous flow equations with an Interaction law based on the 
Hilbert Integral. The latter method has been Incorporated Into the LEWICE code 
so that ice shapes and performance degradation of the airfoil can be predicted. 

The present paper applies the combination of the modified LEWICE code (one 
with the smoothing procedure) and the interactive boundary-layer (IBL) method 
to predict Ice shapes and their effect on airfoil performance. The following 
section presents a brief description of the Interactive boundary-layer proce- 
dure. The results of section 3 allow comparison between calculated and meas- 
ured Ice shapes and the resulting drag increase for a NACA 0012 airfoil. The 
experimental data is due to Olsen, Shaw, and Newton (ref. 9) and were obtained 
in the NASA Icing Research Tunnel (IRT). The measurements were made over a 
large range of conditions which include airspeed and temperature, the droplet 
size and liquid water content of the cloud, and the angle of attack of the air- 
foil. The paper ends with a brief discussion of the Implications of the 
results and a summary of the more important conclusions. 


NOMENCLATURE 

A damping-length constant 

Cp drag coefficient 

Cij Interaction-coefficient matrix 

k s equivalent sand-grain roughness 
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k + dimensionless sand-grain roughness 


L mixing length 

T s static air temperature 

u external velocity 

e 

u° Invlscld velocity 

e 

u friction velocity 

X 

u K Invlscld velocity which contains the displacement effect from a previous 

e sweep 

V n blowing velocity 

Vo, airspeed 

x distance along the surface 

y distance normal to the surface 

y + a Reynolds number, y u^/V 

Su e perturbation velocity due to viscous effects 

S* displacement thickness 

k universal constant, also used as a sweep parameter 
v kinematic viscosity 


2.0 DESCRIPTION OF THE INTERACTIVE BOUNDARY-LAYER METHOD 

The original LEWICE code Is described In detail In reference 4 and an 
operator-free procedure Incorporated Into this code In order to avoid the 
occurrence of the multiple stagnation points caused by the formation of Irregu- 
lar surfaces on the Ice shape Is discussed in reference 6. For this reason, 
the following description of the interactive boundary-layer procedure is brief. 
It couples the solutions of invlscld and viscous flow equations so as to ensure 
that each Influences the other. The Invlscid flow equations are solved by a 
panel method in which the airfoil and ice shapes are defined by a set of points 
where neighboring points are connected by straight-line panels which each have 
source density and vorticity. The vortlcity strength of each panel is the same 
so that vortlcity is defined by a total strength, adjusted to satisfy the Kutta 
condition. The source strengths have independent values on each panel and 
these are adjusted, by solving a set of simultaneous linear equations, to sat- 
isfy the normal-velocity boundary condition at the midpoints of the panels. In 
the strictly Invlscld case this condition requires that the total normal veloc- 
ity, freestream plus body sources and vortices, should vanish. When the bound- 
ary layer Is simulated, the desired normal velocity, v n . Is finite and equals 
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the derivative along the surface of the product of tangential velocity and dis- 
placement thickness, d/dx(u e $*). It Is known that this surface blowing distri- 
bution displaces the dividing streamline outward from the surface of the air- 
foil to the location of the displacement thickness. Experience has shown that 
best results are obtained when surface pressures are calculated and the Kutta 
condition applied on the displacement surface rather than on the surface 
panels. 


The boundary-layer equations for two-dimensional external steady Incom- 
pressible flows are well known and are solved with the Reynolds shear stress 
term modeled with the Cebecl-Smith eddy-viscosity formulation (ref. 10). For 
the external velocity distribution specified by the panel method, u°(x), and 

with $u x representing the perturbation velocity due to viscous effects, the 
edge boundary condition is written as 

u (x) = u?(x) + Su (x) (1 ) 

e e e 


where, for the interaction region confined to the range x a < x < x^, which is 
often taken to include the airfoil chord length plus two chord lengths from the 
tral 1 ing edge , 


Su (x) 
e 



( 2 ) 


with d(u e S*)/do corresponding to the blowing velocity. 

In this form, equation (1) provides an outer boundary condition for the 
viscous-flow calculations which represents the viscous/inviscld interaction and 
can be generalized to the form 


U ( X ) = u*(x) 
e e 




(3) 


where u*(x) corresponds to the inviscid velocity distribution which contains 

the displacement thickness effect, (S*)*, computed from a previous sweep. Here 
denotes the interaction-coefficient matrix, which is obtained from a dis- 
crete approximation to the Hilbert integral. 


Ice on airfoils can introduce substantial geometric changes to their lead- 
ing edges in a short period and cause rapid variations in the flow properties. 
As a result, the inviscid and viscous flow calculations may have difficulty in 
producing satisfactory solutions. For the boundary-layer calculations, the 
Iced airfoil Is regarded as a smooth or a rough surface obtained by covering 
the leading-edge region with a "blanket" as shown in figure 2. It also makes 
use of a continuation method in which the initial calculations are performed 
for the smooth airfoil and subsequent ones for a series of shapes that fall 
between the "smooth" and iced airfoils; For each shape, the blowing velocity 
is computed from 
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v n - 3 T [u (6* - 6) 
n dx L e J 


( 4 ) 


where S* corresponds to the displacement thickness obtained from the 
boundary-layer solutions for the shape whose geometrical difference from the 
"smooth" airfoil Is 6^^(x) and where the S*-surface Is outside the singular- 
ity surface. This allows the viscous effects to be Incorporated Into the 
1 nvl sc 1 d flow solutions gradually, at each time step, thus reducing the sensi- 
tivity of the viscous flow solutions to the rapid changes in the pressure dis- 
tribution near the leading edge. For further details, see reference 6. 

The numerical solutions of the boundary-layer equations, written In trans- 
formed variables, are obtained with the box method for both standard (pre- 
scribed pressure distribution) and Interactive methods. This second-order 
finite-difference method has been used extensively by Cebecl and his associates 
for a wide range of flows (ref. 10). An Inverse form of the equations Is used 
to obtain the solutions with separation and the FLARE approximation, In which 
the convective term u(9u/3x) Is set equal to zero In the recirculating region, 
Is employed. The nonlinear system of algebraic equations which results from 
the finite-difference approximations Is linearized by Newton's method and 
solved by a block elimination procedure (ref. 10). It should be. noted that the 
mixing length expression of the Cebecl-Smlth model (ref. 8) has been modified 
to deal with surface roughness such as that associated with Ice. This was done 
by modifying the mixing length and wall-damping expression, that Is 

L = <(y + Ay) { 1 - exp[-(y + Ay) /A ] } (5) 

where Ay Is a function of an equivalent sand-grain roughness k s . In terms of 
dimensionless quantities, with k* = k*u^/v and Ay + = AyuWv, 



5 < k+ < 70 


70 < k* < 2000 


( 6 ) 


The equivalent sand-grain roughness for Ice Is determined from the expressions 
used In the LEWICE code, as discussed In the following section. 


3.0 RESULTS AND DISCUSSION 

The computer program described In reference 6 has the option of computing 
the flowfleld without and with viscous effects. Studies conducted with this 
code show that the viscous effects do not have a pronounced effect on the pre- 
diction of Ice shapes when the flowfleld Is computed with the Interactive 
boundary-layer procedure rather than with the Invlscld method alone. Even 
though the viscous flowfleld differs significantly from that computed with a 
panel method, especially near the leading edge, the computed shapes remain 
relatively unchanged because of the Insensitivity of the current heat transfer 
model to viscous corrections. As a result, all Ice shape calculations, pre- 
sented herein, are performed with the Invlscld panel method. 
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Before we compare the calculated and measured results. It Is useful to 
review the extensive data of Olsen et al . (ref. 9) for a 0.53-m chord NACA 0012 
airfoil. The experimental data Includes Ice shapes and the resulting drag 
coefficients and were obtained over a range of air temperature, airspeed, air- 
foil angle of attack, spray time, liquid water content and droplet size. The 
data are very Informative and helpful not only In understanding the Ice struc- 
ture and the way It forms, but It Is also very useful In the development of 
computer codes such as the LEWI CE and IBL codes to predict the Ice accretion 
and resulting drag. 

The experiments of Olsen et al . encompass a temperature range that gives rise 
to both rime and glaze Ice. Figure 3(a) shows the repeatability of the Ice 
shape and resulting drag coefficient for a typical rime Ice shape accreted at 
-26 °C. Olsen et al. reported that the Ice shapes and resulting drag coeffi- 
cients repeated quite well and the scatter In the drag coefficient was compa- 
rable to that observed with the clean airfoil data shown In figure 4. 

Figure 3(b), corresponds to a different Icing condition and the rime Ice shape 
and draj coefficient, Cp, again repeated well. The variations for both sets of 
Op data were about ±5 percent, which was close to the variation noted for the 
clean alrfol 1 data. 

Similar comparisons are shown In figures 5(a) and (b) for two glaze Ice 
shapes reported In reference 9. As can be seen, there Is a much larger varia- 
tion In the Ice shapes and drag coefficient for two cases on figure 5(a) and 
one case on figure 5(b). The variations for both sets of Cq data are about 
±1.5 percent of the average values of Cp which Is larger than those for the 
rime shapes or for the clean airfoil. Poor repeatability has also been noted 
In other airfoil tests with glaze Ice. Olsen et al . stated that "no certain 
explanation for the poor repeatability of glaze Ice shapes Is available at this 
time." 


The heat transfer model used in the LEWICE code makes use of an equivalent 
sand-grain roughness, k s , expressed as a function of liquid water content 
(LWC) , static air temperature (T s ), and airspeed (Vc) In order to determine the 
Ice shapes. With c denoting the airfoil chord and <k s / c ) base = 0.001177, it 
Is expressed In the following form 


k $ /c 


i — - 

U 

c/1 

1 


k s /c 

(-) 

v 'case 

m 

[ (k s /c) baseJ 

« 

LWC 

, <k s /c> base. 

T, 

[<V c) basJ 


(7) 


where each sand-grain roughness parameter Is given by 


V c 


<V c Wj 


0.5714 + 0.2457 (LWC) + 1.257KLWC)' 


(8a) 


LWC 


k s /c 


<k s /c) baseJ 


0.047T $ - 11.27 


(8b) 
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(8c) 


y. = 0.4286 + 0.00441 39V 

(K-/C). 00 

L s basejy 
00 

These expressions are empirical, are based on experimental data reported In 
reference 4 and do not account for the effect of time on the Ice roughness. 

The experimental data of Olsen et al . shows that, for the glaze Ice condition, 
the roughness Increases with time, rapidly at first then more slowly. Their 
data also shows that rime Ice Is never as rough as glaze Ice. 

The experimental data of Olsen et al . was obtained at two airspeeds, V*,, 
corresponding to 58 and 94 m/s. The data shows that the roughness Is nearly 
independent of V® of 58 and 94 m/s and that the equivalent sand-grain rough- 
ness parameter k s Is also a function of the median volume droplet (MVD) size, 
as well as a function of the parameters In equation (7). Based on these exper- 
imental observations, It Is plausible to write equation (7) as 



Some numerical studies conducted with the LEWICE code showed that the calcu- 
lated results agreed better with experiment If the roughness parameter for 
was taken as a constant 0.6839 corresponding to V® of 58 m/s. Therefore, all 
subsequent calculations were performed with this assumption. 

An appropriate sand-grain roughness parameter k s Is also required In the 
turbulence model In the IBL code In order to compute the boundary-layer devel- 
opment on the airfoil and In the wake. It Is plausible to conjecture that this 
parameter, as In the heat transfer model, Is a function of LWC, V®, T s and 
MVD. For this reason, calculations with k s corresponding to that used for 
the heat transfer model, equation (9) were performed with the IBL code for an 
Ice shape determined with the LEWICE code. As expected, the drag coefficients 
calculated from the wake velocity profiles Indicated that the results were 
sensitive to the magnitude of the roughness parameter. Several calculations 
with the value of k s of equation (9) multiplied by a constant equal to 2 
yielded the best agreement with experiment. As a result, all boundary-layer 
calculations In the IBL code used a roughness parameter, <k s ) i bl * Qlven by 

< k s> IBL - 2<*<s ) Eq.<9) (11) 
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The following subsection presents the results for the Ice shapes determined 
with the LEWICE code and obtained for a range of conditions which Include V®, 
T s , MVD and LWC, all for a given angle of attach of 4°, and compares them 
with the measured shapes. As stated earlier, these calculations were performed 
with the flowfield computed with the inviscid panel code. Comparisons between 
the calculated and measured drag coefficients, again for a range of conditions 
as above plus the angle-of-attack effect, are presented in subsection 3.2. In 
this case, the flowfield calculations on the airfoil and in the wake were per- 
formed In an Interactive manner for the given calculated ice shapes, as dis- 
cussed In section 2.0. 


3.1 Comparison Between Calculated and Measured Ice Shapes 

Calculations were performed with the LEWICE code in order to determine the 
effects of air temperature, T s , liquid water content LWC and droplet size on 
ice shapes for the experimental data of Olsen et al . (ref. 9) The calculated 
and experimental results are shown by solid and dotted lines, respectively, In 
figures 6 to 12. Those shown in figures 6 and 7 Indicate the effect of air 
temperature on the Ice accretion and. In general, show promising agreement with 
data for two airspeeds. At low temperatures, when the Ice shapes correspond to 
rime Ice, the comparison between calculated and measured ice shapes are good. 
Except for one case in each figure, the calculations were performed without 
numerical difficulties for the stated times of ice accretion In the experiments 
which were 8 min for those in figure 6 and 6.2 min was for those In figure 7. 

The results in figures 8 and 9 show the effect of liquid water content on 
Ice shapes. Those In (a) and (b) of each figure are for values of LWC of 1.0 
and 1.3, and are In good agreement with experiment although the ice shapes of 
(b) are not as good as those of (a). For higher values of LWC, figures 8(c) 
and 9(c), the agreement between measured and computed Ice shapes worsens. The 
results are not unexpected because the correlations for the roughness parameter 
for LWC used In reference 4 are for a range of LWC up to 1.0 g/m^. Further- 
more, the experimental data of Olsen et al. (see also fig. 10, taken from 
ref. 9) Indicate that, for values of LWC ranging from 1.0 to 2, the effect of 
LWC on the ice shape Is negligible and this suggests that the roughness parame- 
ter for LWC should be nearly constant. The calculations for the highest value 
of LWC also Indicate some numerical difficulties with solutions breaking down 
after Ice accretion time of 6 min, rather than the specified time of 8 min. A 
better correlation may improve the predicted Ice shapes and avoid numerical 
difficulties. 

Figures 11 and 12 show the effect of droplet size on the Ice shapes, 
again for two airspeeds. These results are generally in good agreement with 
data and are much better than those which did not Include equation (10) in the 
cal culations. 


3.2 Drag Coefficients 

At first the calculations were performed to investigate the effects of 
droplet size, liquid water content, and air temperature on the total drag coef- 
ficients of the airfoil. The ice shape determined at an angle of attach of 4° 
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with the procedure described above was fixed In the calculations. Tables I 
to III show the results obtained In this manner. 

Tables 1(a) and (b) show the effect of air temperature on drag coefficient 
for Ice accretion times of 8 and 6.2 min for two airspeeds, respectively, 
together with the variation of the equivalent sand-grain roughness parameter 
k s with temperature. Figure 13 shows the variation of the drag coefficient as 
a function of total air temperature and corresponds to the results of table I. 
As can be seen, at lower temperatures where the Ice accretion leads to the for- 
mation of rime Ice, the computed drag coefficients are lower than the measured 
ones but have the same constant level (around 0.013) as the experimental values 
(around 0.020). At higher temperatures, the calculated drag coefficient shows 
a dramatic increase followed by a sudden decrease. In general, the drag coef- 
ficients of those iced airfoils, which are of the glaze type, are in good 
agreement with experiment despite the breakdown of the IBL calculations at 
T s = 2.6655 K in table 1(b). In this case, the calculations capture the 
Increase In the drag coefficient but not the peak value which may correspond 
to stall or post-stall conditions. 

Tables 11(a) and (b) show the effect of liquid water content on drag coef- 
ficient. Except at the high values of LWC, the calculations are performed 
without numerical difficulties and lead to satisfactory results with experimen- 
tal data. Perhaps, as indicated in subsection 3.1, the calculated drag coeffi- 
cients can be Improved if the predicted ice shapes are brought into a better 
agreement with measured ones by improving the roughness correlations for LWC. 

Tables 1 1 1 ( a) and (b) show the effect of droplet size on drag coefficient, 
again for two airspeeds. The results in table I I I (a) .contain relatively lower 
drag coefficients indicating smaller viscous effects. For the three values of 
the droplet size, the computed values are in reasonable agreement with experi- 
ment. The experimental drag coefficients of table 111(b), on the other hand, 
are very large for values of MVD > 20, suggesting the existence of increased 
flow separation and stall and/or post-stall conditions. For these conditions, 
the computed values are in poor agreement with data. 

Figures 14 and 15 show the variation of the drag coefficient of the air- 
foil with angle of attack for a given ice shape determined at a = 4°. The 
results in figure 14 are for rime ice and indicate remarkably good agreement 
between calculations and experiment. The Increase In the drag coefficient due 
to the Ice is well represented and the turbulence model with equivalent sand- 
grain roughness parameters as given by equation 10 allows the calculations to 
follow the experimental trend. 

The results In figure 15 correspond to glaze ice and are again in remark- 
ably good agreement with experimental data except at higher angles of attack 
where the code broke down due to the occurrence of stall. It should be noted 
that the calculations were performed for the ice shape with LWC = 1.3 g/nw 
rather than 2.1 g/m^ since, according to figure 10, the Ice shapes are not sen- 
sitive to LWC. In figures 8 and 9, it has already been shown that predicted 
Ice shapes are not in good agreement with data at the higher values of LWC and 
Indicates the need for Improved roughness correlations. Since the accuracy of 
the drag coefficient is a strong function of the Ice shape, it is necessary to 
use a shape that fits the experimental data as well as possible. 
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4.0 CONCLUDING REMARKS 


The extensive results of the previous section show that the interactive 
boundary- layer method of section 2.0 and the modified LEWICE code provide a 
satisfactory basis for the calculation of ice shapes and their effect on air- 
foil performance. The results encompass a wide range of icing conditions and 
angles of attack of direct relevance to engineering practice so that the method 
provides a basis for extension to the representation of icing on wings, other 
lifting bodies and engine Intakes. 

Improvements can be made to the present approach, for example to the cor- 
relation for equivalent sand-grain roughness and to the representation of the 
drag coefficient when the ice accretion gives rise to stall or post-stall con- 
ditions. In the former case, additional measurements and calculations are 
required and in the latter a combination of local flow measurements and calcu- 
lations to examine the sensitivity of the aerodynamic flow characteristics to 
the time-change in ice shapes. The improvements should be pursued but should 
not preclude the immediate extension of the calculation method to deal with 
three-dimensional lifting bodies and the development of a better heat transfer 
formulation so that the ice accretion model can also include the effects of 
surface heating on ice formation. 
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TABLE 1(a). - EFFECT OF AIR TEMPERATURE ON DRAG COEFFICIENT 
[Droplet size, MVD, 20 pm.] 

(a) Airspeed, V x , 58 m/s; liquid water content, 


LWC, 1.3 g/m^; ice accretion time, 8 min 


Static air 
temperature , 

T s- 

K 

Drag coefficient, Cg 

Equi val ent 
sand-grain 
roughness , 

k s* 

mm 

Cal culated 

Experimental 

245.35 

0.01279 


0.334 

253.35 

.01346 


.815 

259.25 

.01431 

.02072 

1.170 

266.45 

.05075 

.06036 

1.602 

269.25 

.02873 

.02807 

1 .770 

270.35 

.02105 

.02647 

1 .837 


(b) ""Airspeed , 94 m/s; liquid water content^ 


LWC, 1,05 g/m^; ice accretion time, 6.2 min 


Static air 
temperature , 

V 

K 

Drag coefficient, Cg 

Equi val ent 
sand-grain 
roughness , 
k$ , 
mm 

Cal cul ated 

Experimental 

242.65 

0.01143 

0.0238 

0.126 

256.55 


.0370 

.740 

260.95 



. .934 

266.55 

(a) 

.0756 

1.182 


breakdown . 


TABLE II. - EFFECT OF LIQUID WATER CONTENT ON DRAG COEFFICIENT 

[Droplet size, MVD, 20 pm; airspeed, V^, 

58 m/s; ice accretion time, 8 min.] 


(a) Static air temperature, T s , 253.35 K 


Liquid water 
content, 
LWC, 

g/m^ 

Drag coefficient, Cg 

Equivalent 
sand-grai n 
roughness , 
k$ , 
mm 

Cal culated 

Experimental 

1 .0 

0.01398 

0.0212 

mSm 

1.3 

.01346 

.0246 


2.0 

(a) 

,0312 

Bl 

(b) Static air temperature, T $J 263 

.48 K 

Liquid water 

Drag coefficient, Cn 

Equivalent 

rnntpnt’ 



sand-grai n 

LWC, 

Cal culated 

Experimental 

roughness , 

g/m^ 



k s , 




mm 


0,02976 

0.0262 



.04585 

.0307 


■V 

(a) 

.0456 

Hfl 


breakdown. 
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TABLE III. - EFFECT OF DROPLET SIZE ON 
DRAG COEFFICIENT 

(a) Airspeed, V , 58 m/s; liquid water content, 

LWC, 1.3 q/m, static air temperature, T $ , 
245.35 K; ice accretion time, 8 min 


Droplet | Drag coefficient, Cq 

Calculated Experimental 


Equivalent 
sand-grain 
roughness , 


0.01283 

.01279 

.01285 


0.0121 

.0193 

.0196 


Droplet 

size 

Drag coefficient, Cq 

Equivalent 
sand— gra i n 

mvd! 

pjn 

Calculated 

Experimental 

roughness , 
k s , 

mm 

14 

0.01578 

0.03090 

0.934 

20 

.04549 

.09510 

.934 

26 

,05356 

.12660 

.747 


Example of real ice accretion frnm Ref 


BLANKET (SMOOTH AIRFOIL) 




IK 


Boundary-layer model used In the viscous 
flow calculations . 





















C D 

DIFFERENCE 


FROM 


AVERAGE 


per (Pent 

0.01941 

0.9 

, 02009 

4.5 

. 01812 

-6.0 

— .01930 

.4 


% VARIATION, + 4.0 PERCENT 



01053 

3.0 

01077 

5.3 

01003 

-10 

00957 

-6.4 


* VARIATION, ♦ 5.0 PERCENT 



(a) MVD, 20ym: LWC , 1.3g/m^: (b) MVO, 12ym: LWC, 1.08 g/m^: 

time, 8 min, time, 5 min. 

Fig. 3. Repeatability of Ice shape and drag for rime Ice shapes. Total temp- 
erature, -26°C; airspeed, 58 m/s; 0.53 m-chord NACA 0012 airfoil at 4® 
(Ref. 9). 


O DATA FOR 0.53-m-CHORD PRODUCTION 



Fig. 4. Comparison of measured clean airfoil data 
with published data for the NACA 0012 airfoil (Ref. 9). 
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C D 

DIFFERENCE 

FROM 

AVERACE 


PERCENT 

0, 033S2 

•3.4 

.02767 

•21.0 

. 03729 

+6.5 

.04134 

+18.0 


* VARIATION, + 16.0 PERCENT - 



(a) MVD , 20ym: LWC, 2.1g/m 3 : 
time, 5 min. 


. 02960 

-3.0 

1 

1 

1 

\ 

<3 

-4.4 

.02949 

-3.3 

. 03SS4 

+20.7 

% VARIATION, + 14.0 

PERCENT 



(b) MVD, 20ym: LWC, 1.3 g/m 3 : 
time, 8 min. 


Fig. 5. Repeatability of ice shape and drag for glaze Ice shapes. Total temp- 
erature, -8°C; airspeed, 58 m/s; 0.53 m-chord NACA 0012 airfoil at 4° 
angle (Ref. 9). 
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Predicted 

: Measured 



(d) 


Effect of air temperature T s (°C) on Ice shapes for fixed droplet 
size (MVO = 20 pm), air velocity (V a = 94 m/s), liquid water 
content (LWC = 1.05 g/m3): (a) -30.5°, (b) -16.6°, (c) -12.2°, ( 

-6.6°. All calculations are for 6.2 min. except for that of (d) 
which Is for 4.13 min. 



rt- 






Predicted 

Measured 



(c) 


Effect of liquid water content, LWC (g/m 3 ), on Ice shapes for fixed 
airspeed (V ra = 58 m/s), temperature (T s = -9.67°C), droplet size 
(MVO = 20 pm): (a) 1.0, (b) 1.3, (c) 1.6. All calculations are for 

8 min. 



(a) (b) 


Effect of LWC on the Ice shape and section drag. Airspeed, 58 m/s 
MVO, 20 pm; time, 8 min at 4° angle of attack (Ref. 9). 






Predicted 

— Measured 



Fig. 12. Effect of droplet size, MVD (pm), on Ice shapes for fixed airspeed 
(Vo, * 94 m/s), temperature (T s = -12.20°C), liquid water content 
(LWC = 1.05 g/m 3 ): (a) 14, (b) 20, (c) 26. All calculations are 
for 6.2 min. 


0.10 r V=58m/s,lWC=>1.3Cg/m ,t=8min 


• CALCULATED 
O EXPERIMENT 



0.00 1 L - — 1 ' ‘ 1 

-30 -25 -20 -15 -‘O -5 0 

TOTAL TEMPERATURE, iO 


Fig. 13. Variation of the drag coefficient with 
total temperature at two air speeds. (Line is a fit 
for experimental data.) 
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Fig. 14. Variation of the drag coefficient of the 
airfoil with angle of attack with Ice shape deter- 
mined at a = 4° for Vo, = 58 m/s, T s = -27.80°C, LWC 
= 1.0 g/m 3 , MVD = 12 ym, Ice accretion time = 
5 min. 



Fig. 15. Variation of drag coefficient of the air- 
foil with angle of attack with Ice shape corres- 
ponding to glaze Ice determined at a = 4° for 
V„ = 58 m/s, T s = -9 . 67°C , LWC = 1.3 g/m 3 , 

MVD = 20 ym. Ice accretion time = 5 min. 
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